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ABSTRACT 



Nbody simulations are used to examine the consequences of Neptune's out- 
ward migration into the Kuiper Belt, with the simulated endstates being com- 
pared rigorously and quantitatively to the observations. These simulations con- 
firm the findings of Chiang et al. (2003), who showed that Neptune's migration 
into a previously stirred-up Kuiper Belt can account for the Kuiper Belt Objects 
(KBOs) known to librate at Neptune's 5:2 resonance. We also find that capture 
is possible at many other weak, high-order mean motion resonances, such as the 
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11:6, 13:7, 13:6, 9:4, 7:3, 12:5, 8:3, 3:1, 7:2, and the 4:1. The more distant of these 
resonances, such as the 9:4, 7:3, 5:2, and the 3:1, can also capture particles in 
stable, eccentric orbits beyond 50 AU, in the region of phase space conventionally 
known as the Scattered Disk. Indeed, 90% of the simulated particles that persist 
over the age of the Solar System in the so-called Scattered Disk zone never had 
a close encounter with Neptune, but instead were promoted into these eccentric 
orbits by Neptune's resonances during the migration epoch. This indicates that 
the observed Scattered Disk might not be so scattered. This model also produced 
only a handful of Centaurs, all of which originated at Neptune's mean motion 
resonances in the Kuiper Belt. However a noteworthy deficiency of the migration 
model considered here is that it does not account for the observed abundance of 
Main Belt KBOs having inclinations higher than 15°. 

In order to rigorously compare the model endstate with the observed Kuiper 
Belt in a manner that accounts for telescopic selection effects, Monte Carlo meth- 
ods are used to assign sizes and magnitudes to the simulated particles that survive 
over the age of the Solar System. If the model considered here is indeed represen- 
tative of the outer Solar System's early history, then the following conclusions are 
obtained: (i) the observed 3:2 and 2:1 resonant populations are both depleted by 
a factor of ~ 20 relative to model expectations; this depletion is likely due to un- 
modeled effects, possibly perturbations by other large planetesimals, (m.) the size 
distribution of those KBOs inhabiting the 3:2 resonance is significantly shallower 
than the Main Belt's size distribution, {iii.) the total number of KBOs having 
radii R> 50 km and orbiting interior to Neptune's 2:1 resonance is iV ~ 1.7 x 10^; 
these bodies have a total mass of M ~ 0.08(p/l gm/cm^)(p/0.04)~^/^ M® assum- 
ing they have a material density p and an albedo p. We also report estimates of 
the abundances and masses of the Belt's various subpopulations (e.g., the reso- 
nant KBOs, the Main Belt, and the so-called Scattered Disk), and also provide 
upper limits on the abundance of Centaurs and Neptune's Troj well as up- 

per limits on the sizes and abundances of hypothetical KBOs that might inhabit 
the a > 50 AU zone. 

Subject headings: Kuiper Belt — solar system: formation — methods: N-body sim- 
ulations 
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1. Introduction 

The Kuipcr Belt is the vast swarm of small bodies that inhabit the outer Solar System 
beyond the orbit of Neptune. These Kuiper Belt Objects (KBOs) that inhabit this Belt are 
relics of the solar system's primordial planetesimal disk — they are bits of debris that failed to 
coalesce into other large planets. The Kuiper Belt is also of great interest since it preserves 
a record of the outer Solar System's early dynamical history. This is reflected in the KBOs' 
curious distribution of orbits, which suggest that there was considerable readjustment of the 
Solar System's early architecture. The possibility that the orbits of the giant planets may 
have shifted significantly (that is, after the solar nebula gas had already dissipated) was first 
demonstrated by the accretion simulations of Fernandez & Ip (1984); they showed that as 
the growing giant planets gravitationally scatter the residual planetesimal debris, they can 
exchange angular momentum with the debris disk in a way that causes the planets' orbits to 
drift. Malhotra (1993b) later showed that an episode of outwards migration by Neptune by 
at least Aa ~ 5 AU could also account for Pluto's pecuhar orbit, which resides at Neptune's 
3:2 resonance with an eccentricity of e ~ 0.25. In this scenario, Pluto's large eccentricity is a 
consequence of it having been captured by Neptune's advancing 3:2 resonance, which pumped 
Pluto's e up as it shepherded the small planet outwards. Further support for this planet- 
migration scenario is provided by the subsequent discovery of numerous other KBOs also 
inhabiting Neptune's 3:2 resonance with eccentricities similar to model predictions (Malhotra 
1995), as well as by more modern Nbody simulations of the orbital evolution of giant planets 
while they are still embedded in a massive planetesimal disk (Hahn & Malhotra 1999; Gomes 
et al. 2004). 

The purpose of the present work is to use higher-resolution simulations to update this 
conventional model of Neptune's migration into the Kuiper Belt. This model's strengths, 
as well as its weaknesses, will be assessed quantitatively by rigorously comparing the sim- 
ulations' endstates to current observations of the Belt. In the following, we execute two 
simulations that track the orbital evolution of the four migrating giant planets plus 10^ 
massless test particles (the latter representing the KBOs) over the age of the Solar system. 
In one simulation the initial state of the Kuiper Belt is dynamically cold {i.e., the particles 
have initial eccentricities and inclinations of e and sini ~ 0.001), while the second simula- 
tion is of a Kuiper Belt that is initially stirred-up a modest amount (i.e., e and sinz ~ 0.1). 
We then use a Monte Carlo method to assign sizes (and hence magnitudes) to the simulated 
KBOs; this allows us to account for the telescopic biases that tends to select those KBOs that 
inhabit orbits that are more favorable for discovery over those KBOs in less favorable orbits. 
Then, by comparing the resulting model Kuiper Belts with the current observational data, 
we rigorously test the planet-migration scenario as well as obtain a more realistic assessment 
of the abundance of KBOs. This analysis will also provide the relative abundance of the 
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Belt's various subpopulations — the resonant KBOs, the Main Belt Objects, the Scattered 
Disk, plus the Centaurs and Neptune's Trojans. 

The paper is organized as follows. Section 2 describes the so-called 'standard model' 
that is considered here, as well as the numerical methods to be employed. Results from two 
simulations of the Kuiper Belt are reported in Sections 3 and 4, while Section 5 examines 
the Kuiper Belt inclination problem. Section 6 details the Monte Carlo model that is then 

used in Sections 7-10 to assess the relative abundance of the Belt's various subpopulations, 
with a final tally of the Belt's total population given in Section 11. Section 12 comments on 
some important unmodeled effects, and Section 13 summarizes the results. 



2. Simulating planet migration 

The MERCURY6 Nbody integrator (Chambers 1999) is used to track the orbital evo- 
lution of the four giant planets plus numerous massless particles. In our simulations, planet 
migration is implemented by applying an external torque to each planet's orbit so that its 
semimajor axis aj varies as 

«.W=a/,,-A,e-*/^ (1) 

where afj is planet j's final semimajor axis, Aj is the planet's net radial displacement, and 
T is the e-fold timescale for planet migration; this form of planet migration was first used 
in Malhotra (1993b). To implement this in MERCURY6, the integrator is modified so that 
each planet's velocity Vj is incremented by the small velocity kick 

A^;^ = l^^e-*/-^. (2) 

with each timestep At. This additional velocity kick is directed along the planet's velocity 
vector, and results in a torque Tj = rrijajAvj/At being applied to each planet. Since 
Tj = dLj/dt where Lj = the planet's angular momentum, these velocity kicks cause the 
planet's orbit to vary at the rate dj = 2ajTj/Lj = (Aj/T)e~*/'^, which then recovers Eqn. (1) 
when integrated. 

The simulations reported below adopt the current planets' masses and orbits as initial 
conditions, except that their initial semimajor axes are displaced by an amount — Aj so that 
the migration torque ultimately delivers these planets into orbits similar to their present 
ones. The free parameters that describe this migration are the planets' radial displacements 
Aj and the migration timescale r. At present, there is only one strong constraint on the 
Aj's, namely, that Neptune's orbit must expand by A^v — 8.7 AU if resonance trapping is 
to account for the KBOs having eccentricities of < e < 0.32 at Neptune's 3:2 resonance 
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(see Appendix A). Another constraint, on the magnitude of Jupiter's inward migration, can 
be obtained from the orbital distribution of asteroids in the outer asteroid belt. Liou & 
Malhotra (1997) show that the severe depletion of the outer asteroid belt can be explained 
if Jupiter migrated inward by at least 0.2 AU, and Franklin et al (2004) show that the 
orbits of the Hilda asteroids at Jupiter's 3:2 resonance are consistent with Jupiter having 
migrated inwards by about 0.45 AU. The remaining A^'s for Saturn and Uranus are less well- 
constrained, but stability considerations require them to be neither too large nor too small. 
With this in mind, our simulations adopt the following values for the A^'s: Aj = 0.900 AU 
for Jupiter, A5 = 2.09 AU for Saturn, A[/ = 4.52 AU for Uranus, and Aat = 8.70 AU for 
Neptune. All of the simulations reported here also employ a planet-migration timescale of 
T = 10^ years. This value is supported by the self-consistent Nbody simulations by Hahn & 
Malhotra (1999) of the giant planets' migration while they are embedded in a planetesimal 
disk. Those simulations show that a planetesimal disk having a mass M/) ~ 50 spread 
over 10 < a < 50 AU will cause Neptune's orbit to expand Aa ~ 7 AU over a characteristic 
timescale of r ~ 10^ years (see also Gomes et al. 2004). 

Wc note that the orbital evolution adopted here is constructed so that the migrating 
planets' eccentricities are always comparable to their present values, and that the migration 
proceeds along nearly circular orbits. But this particular choice for the planets' eccentricities 
is merely a simplifying assumption since we do not know the e-evolution of the giant planets 
during the migration epoch. For instance, it is possible that dynamical friction with the 
particle disk would have conspired to keep the planets' eccentricities low, but there may also 
have been other transient protoplanets roaming about the outer Solar System, and their 
perturbations would tend to pump up the planets' eccentricities. Given the uncertainty in 
the relative rates of these effects, wc adopt the simplest possible model, one that assumes 
that the giant planets' eccentricities were always comparable to their present one. However 
alternate migration schemes are possible; for instance, Tsiganis et al. (2005) consider a 
scenario where the giant planets pump up their eccentricities as they pass through mutual 
resonances. But it is uncertain as to whether this possible history would have altered the 
bulk properties of the Kuiper Belt, and it is not considered here. 

In order to enforce migration in nearly circular orbits, our simulations have Jupiter mi- 
grating outwards Aj = 0.9 AU, whereas other self-consistent simulations show that Jupiter 
usually migrates inwards a small amount (Hahn & Malhotra 1999). Note that this choice 
avoids having Jupiter approach the 5:2 resonance with Saturn, which tends to excite the 
planets' eccentricities above current levels. But that eccentricity excitation might then have 
been damped back to current levels by dynamical friction with the particle disk, but that is 
a phenomena that goes unmodeled in our massless particle disk. We simply avoid this event 
altogether by instead having Jupiter migrate outwards a modest amount. But this not a 
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concern here since our interest is in the Kuiper Belt, whose endstate is not hkely to preserve 
any memory of whether Jupiter migrated shghtly inwards or outwards. The remaining A^'s 
are similarly chosen to avoid all major resonances, and the simulated planets' semimajor 
axes aj{t) are also shown in Fig. 1. 

In the simulations described below, the model Kuiper Belt is initially composed of 10^ 
massless particles having semimajor axes randomly distributed over 20 < a < 80 AU with 
a surface number density that varies as a~^. The inner edge of this particle disk is 1.4 AU 
inwards of Neptune's initial semimajor axis, and the disk extends well beyond the outer 
reaches of the observed Main Belt. All simulations described here use a timcstcp of At = 0.5 
years, which is sufficiently short to accurately evolve particles in eccentric orbits down to 
perihelia as low as Qmin — (10 At /I yr)^/^ AU ~ 3 AU without suffering the perihehon 
instability^ described in Ranch & Holman (1999). Of course, particles can still achieve 
orbits having perihelia lower than qmin, but their orbits will not be calculated correctly in 
our simulations. However this is of little consequence since these planet-crossing bodies have 
very short dynamical lifetimes and are quickly removed from the system anyway. 



3. Migration into a dynamically cold Kuiper Belt 

Accretion models have shown that the observed KBO population must have formed in an 
environment that was initially dynamically cold, that is, the known KBOs must have formed 
from seeds that were in nearly circular and coplanar orbits with initial e's and sin f's that were 
< 10~^ (Stern 1996; Kenyon & Luu 1999). In anticipation of this, many investigations of the 
dynamical history of the Kuiper Belt have adopted initial KBO orbits that are dynamically 
cold (e.g., Malhotra 1993b, 1995; Duncan et al. 1995; Yu & Tremaine 1999; Chiang & Jordan 
2002). 

Figure 2 shows the results of a simulation of Neptune's smooth migration into a dy- 
namically cold swarm of massless Kuiper Belt objects having initial e's that are Raylcigh 
distributed about a mean value (e) = 0.001, and initial inclinations similarly distributed 
with a mean (sini) = (e)/2. This system is evolved for t = 5 x 10^ years. As is well known 
from previous studies, Neptune's smooth migration is very efficient at inserting particles into 
the planet's mean-motion resonances, principally the 2:1, 5:3, and the 3:2. However, it is 
also well recognized that the endstate of this ideahzed model differs from the observed KBO 
orbits (the red dots in Fig. 2) in several ways. For example, one prominent discrepancy 
is that the 2:1 resonance is densely populated with simulated particles while only sparsely 



This difficulty is overcome by the algorithm of Levison & Duncan (2000). 
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populated by observed KBOs. 

However the discrepancy that is most important to this discussion hes in the 44 < a < 47 
AU zone between the 7:4 and the 2:1 resonances, which is the outer half of the Main Belt 
that is conventionally defined as the region between the 3:2 and 2:1 resonances. Although 
these simulated particles managed to slip through the advancing 2:1 resonance, they still 
reside in orbits that arc only modestly disturbed with e ~ 0.05 and i ~ 0.5°, whereas the 
observed KBOs inhabit orbits that are considerably more excited. Thus Neptune's smooth 
migration into a dynamically cold Kuiper Belt is unable to account for the Belt's stirred-up 
state. 

Evidently, some other event has also disturbed the Kuiper Belt, and this stirring event 
may have taken place prior to or after the onset of Neptune's migration. However Section 
4 provides reason to believe that this stirring event occurred before the onset of Neptune's 
migration into the Kuiper Belt. 

4. Migration into a stirred— up Belt 

To examine the effects of Neptune's migration and its resonance sweeping of a previously 
stirred-up Kuiper Belt population, we repeat the numerical integrations with 10^ simulated 
KBOs, but with initial e's Rayleigh distributed about a mean value of (e) = 0.1 and initial i's 
distributed similarly about a mean value of (sini) = (e)/2. However this time the simulation 
is evolved for the age of the Solar System, 4.5 Gyrs, with Fig. 3 showing the resulting Kuiper 
Belt endstate. 

First, we note that in this case we find an outer Solar System that is far more depleted in 
transient particles like Centaurs (which are scattered particles having semimajor axes interior 
to Neptune) and Scattered Disk Objects [which are particles that were lofted into eccentric 
Neptune-crossing orbits due to a close-encounter with Neptune (Duncan & Levison 1997)]; 
those bodies usually reside in orbits having perihelia q between the q — 30 and q — AO 
AU curves seen in Figs. 2 — 3. This difference is primarily due to the simulation's longer 
integration time. 

Another prominent difference with the 'cold belt' simulation is that Neptune's weaker 
higher-order resonances, such as the 3:1 and 5:2, are considerably more efficient at captur- 
ing particles when Neptune migrates into a hot disk, a phenomenon that was first noted by 
Chiang et al. (2003). This result was rather surprising, because low-order resonance capture 
theory theory predicts a generally lower capture probability for particles having higher ec- 
centricities (Borderies & Goldreich 1984; Malhotra 1993a). However, a careful examination 
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of the theory of adiabatic resonance capture (e.g., Dermott et al. 1988) shows that there are 
two reasons for this result. (1) The higher order resonances have capture probabihties that 
drop off more slowly with eccentricity than first order resonances: although the 1st order 
resonance capture probability varies as ~ e~^/^, the 2nd order resonance capture probability 
falls off more slowly as ~ while the 3rd order resonance capture probability varies as 
~ e~^/^. (2) The threshold migration speeds for adiabatic resonance capture are also lower 
for the higher order resonances, and they also depend more strongly upon the initial eccen- 
tricities. For capture at a j+k : j resonance, the requirement for adiabatic resonance capture 
is that Neptune's migrate rate, which is d^^ ~ Ajv/r ~ 10~^ AU/yr in these simulations, be 
sufficiently slow, namely, that 

I • I ^ o -r^ k^N a-N /o\ 

\aN\ < SjCjkC ji—^- (3) 

where rriN and Pn are Neptune's mass and orbital period, and Cjk is a function of Laplace 
coefficients. For example, the 5:2 resonance has C23 — 3.3, so the migration speed threshold 
that permits adiabatic resonance sweeping is |djv| ^ 5 x 10^^ AU/yr among particles having 
e ~ 0.1, while the threshold is reduced to |d7v| ^ 5 x 10^^^ AU/yr among particles having 
e ~ 0.001. It is clear then that a dynamically cold particle swarm has no chance of adiabatic 
capture at Neptune's high-order mean-motion resonances, while particles that are stirred 
up to e ~ 0.1 are at least near the threshold for adiabatic resonance capture. And as Chiang 
et al. (2003) point out, the fact that seven eccentric KBOs are known to hbrate at Neptune's 
5:2 resonance also lends support to the pre-stirred Kuiper Belt scenario. 

Another advantage of this stirred-up Kuiper Belt scenario is that it recovers eccentrici- 
ties that are observed to be as large as e ~ 0.2 in the Main Belt that hes between the 3:2 and 
2:1 resonances at 40 < a < 48 AU (the red dots in Fig. 3). This is a feature that the cold 
Belt scenario (Fig. 2) does not account for. Of course. Fig. 3 also shows that the simulated 
Main Belt is densely populated by low-eccentricity particles having e ~ 0.05 at a ~ 47 AU, 
whereas the observed Kuiper Belt is only sparely populated here. But Section 7 will show 
that there are a variety of possible explanations for this discrepancy — such a change in the 
KBO size distribution, or perhaps an outer edge in the primordial Kuiper Belt. 

Figure 3 also shows that trapping at the distant high-order resonances like the 5:2 and 
3:1 is quite effective at promoting bodies into eccentric orbits having a > 50 AU and perihelia 
30 ^ ^ 40 AU. This domain is usually regarded as the Scattered Disk. This result then 
suggests the possibility that some of the observed KBOs in the 30 < g < 40 AU zone may 
actually be resonantly trapped bodies that are masquerading as members of the Scattered 
Disk. Of course, particles scattered by Neptune also tend to spend a large fraction of their 
time near resonances due to the resonance sticking phenomenon (e.g., Duncan & Levison 
1997; Malyshkin & Tremaine 1999). Therefore, the discrimination between scattered and 
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resonantly trapped particles must be done carefully. Towards this end, we examine the 
orbital histories of all surviving particles in the shaded zone in Fig. 4 that have e > 0.25 and 
a > 50 AU. The usual definition of being 'in' a j + A; : j mean motion resonance is that the 
particle's resonance angle (pj^, Eqn. (Al), librates about some fixed value with some modest 
amplitude |A0jfc| that is usually < 90°, while the resonance angle for a scattered body that 
is temporarily 'stuck' in a resonance will have a (pjk that circulates over ±180°. However, we 
find that |A0j7;| is not the best discriminant for identifying trapped and scattered particles 
because a small but significant fraction of particles do get trapped at a resonance with a (pjk 
that is either circulating or else librating with a very large amplitude. For some trapped 
particles, this distinction is unclear due to this simulation's infrequent time-sampling that 
occurs every At — 10^ years. 

Rather, a more reliable discriminant between trapped and scattered particles is based on 
Brouwer's integral B, Eqn. (A6). This integral is conserved by resonantly trapped particles 
but is not conserved by scattered particles that are temporarily exhibiting the 'resonance 
sticking' phenomenon. Of the 134 particles that inhabit the gray zone in Fig. 4, only 12, or 
about 10% of these particles, are truly scattered particles whose orbits (a, e) evolve stochas- 
tically about the 30 < g < 40 AU zone; these scattered particles are indicated by the crosses 
in Fig. 4. The remaining particles are resonantly trapped particles, most^ of which preserved 
their B integral to within ~ 3%. 

The orbits of all particles having perihelia g < 40 AU have also been inspected, and 
those resonances inhabited by trapped particles having libration amplitudes \(f)jk\ < 90° are 
indicated by the vertical dashes in the Figure. We find that particles get trapped at a number 
of exotic resonances like the 11:6, 13:7, 13:6, 9:4, 12:5, 8:3, and the 11:4. 

5. Kuiper Belt inclinations 

Inspection of the inclinations shown in Fig. 3 suggests that the smooth migration model 
docs not produce sufficient numbers of bodies in high-inclination orbits. This has been 
recognized in previous studies (Malhotra 1995; Gomes 1997), but has not been quantified. 
However one should not directly compare the simulation's i's to the observed KBO incli- 
nations, since the latter is heavily biased by telescopic selection effects. Note that most 
telescopic surveys of the Kuiper Belt observe near the ecliptic, which favors the discovery 



^However the e's and i's of some resonantly trapped particles will still oscillate with constant a in a 
manner that preserves their .Jacobi integral; this evolution usually occurs after migration has ceased, and 
these particular motions do not preserve B. 
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of lower-i KBOs that spend a larger fraction of their time at lower latitudes (Jewitt & Luu 
1995). To mitigate this selection effect, one should instead consider the ecliptic inclination 
distribution, which is the inclination distribution of objects having latitudes j3 very near the 
ecliptic (e.g.. Brown 2001). The ecliptic inclination distribution for multi-opposition KBOs 
having perihelia g < 42 AU and latitudes \j3\ < 1° (red curve) is shown in Fig. 5, as well 
as the simulated ecliptic inclination distribution for particles from Fig. 3 that are selected 
similarly (black curve). 

The agreement seen in Fig. 5 among bodies having inclinations of < i < 15° shows 
that the smooth migration model can readily recover the Kuipcr Belt's lower inclination 
members. Of course, this agreement is partly due to the particles' initial inclinations being 
distributed around {i) ~ 3°. But Fig. 5 also shows this model to be quite deficient in 
producing sufficient numbers of the high-i bodies having i > 15°. Similar results are also 
obtained among bodies orbiting at Neptune's 3:2 resonance. This is a serious discrepancy, 
since Brown (2001) has shown that there are two inclination-populations in the Kuiper 
Belt: a minor population of low-i having characteristic inclinations of i ~ 3°, and a high-i 
population having i ~ 15° containing about three-quarters of all KBOs. Note that these 
high-i bodies are very underrepresented in Fig. 5 due to telescopic selection effects. 

Of course, Neptune-scattered particles routinely achieve high inclinations of i > 15°; for 
instance, many of the high-i particles seen in Fig. 2 were scattered by Neptune. Could the 
Scattered Disk be a source of the high-i KBOs that are found elsewhere in the Belt? Recent 
Nbody simulations by Gomes (2003) show that a small fraction of these Ncptunc-scattcrcd 
particles can evolve from very eccentric, Neptune-crossing orbits into less eccentric orbits in 
the Main Belt. In Gomes' simulations, this occurs principally at secular and mean motion 
resonances that drive large oscillations in a scattered particle's eccentricity. When a scattered 
particle visits a resonance, it can have its e temporarily lowered and its q raised. If this 
occurs during the planet-migration epoch, this process becomes irreversible and can strand 
Scattered particles in the Main Belt with their high inclinations. Such bodies are identified 
by Gomes as 'evaders' since they are Neptune-crossing bodies that ultimately manage to 
evade Neptune when deposited in the Main Belt. Note, however, that the efficiency of this 
process is quite low, affecting only e ~ 0.002 of all Scattered particles in the simulation that 
is evolved over the age of the Solar System by Gomes (2003). However all of the high-i 
particles seen in our simulation (Fig. 3) achieved their inclinations while temporarily or 
permanently trapped in Neptune's advancing resonances. There were no Neptune-scattered 
evaders having i > 10° that survived in our simulations. 

Despite the evader mechanism's inefficiency, a model can still be constructed that yields 
a KBO inclination distribution that is quite similar to the observed one. For instance, this 
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can be achieved by making the number density of small bodies initially orbiting interior 
to ~ 27 AU about 60 times higher than the density of bodies initially orbiting beyond 27 
AU. As Gomes (2003) show, Neptune's migration through this densely-populated inner disk 
creates the Kuiper Belt's high-? evaders, while the sparse outer disk provides the Belt's low-i 
component. Although this scheme yields an i-distribution that does indeed agree with the 
observations, that success is achieved via a special configuration of the initial particle disk. 

However, Levison & Morbidelli (2003) avoid this problem of special initial conditions by 
assuming that the initial planetesimal disk simply ended at ~ 35 AU. This is the 'push-out' 
model, which argues that most of the Kuiper Belt is a consequence of Neptune's advancing 
2:1 resonance, which can drag bodies outwards to litter the Main Belt with low-i KBOs. 
The Belt's high-i component is then presumed to be due to the evader mechanism. While 
the push-out model remains quite intriguing, it would be interesting to see this scenario 
subjected to greater scrutiny to see whether it can indeed reproduce the Kuiper Belt's 
curious mix of high and low inclination KBOs in a self-consistent manner. 



6. A Kuiper Belt census: comparison with observations 

Figure 6 plots the relative abundance, over time, of the simulated Belt's various dy- 
namical classes among particles having perihelia g < 45 AU. These curves are normalized 
such that the final abundance of the Main Belt (MB), where 40.1 < a < 47.2 AU, is unity. 
Note that this model predicts a 2:1 resonance that is 1.4 times more abundant than the 
Main Belt, and 2.5 times more abundant than the 3:2, while the observations (Fig. 3) show 
a 2:1 that is only sparsely populated. Of course, when comparing the simulated population 
to the observed population, one must first deal with the observational selection effect that 
strongly favors the discovery of larger and/or nearer KBOs. However it is shown below that 
the effects of this bias can be accounted for by using a Monte Carlo method that assigns 
random sizes to the simulated population. This then allows one to make a fair comparison 
of the relative abundances of the simulated and observed populations. 

Begin by letting N{R) = the number of bodies in the simulated population having radii 
exceeding a radius R. Also let ^ be a random number that is uniformly distributed between 
zero and one, and interpret this number as the probability of selecting a body with a radius 
that is smaller than R. This is also equal to the probabihty of not selecting a body of radius 
> i?, so ^ = 1 — N [R] / Ntotai where Ntotai is the total number of bodies in this population. 
Since most small-body populations have a cumulative size distribution N{R) that varies as 
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a power law, adopt 



N{R) = N,^tai{^) (4) 



where Rmin is the radius of the smallest member of the swarm. Then R/ R^^n = (1 ~0 ^^"^^ 
but 1 — ^ can be replaced with ^ since these random numbers have the same distribution, so 

^(0 = (5) 

Equation (5) is then used to generate random sizes for the simulated population of Fig. 3 
that have apparent R-band magnitudes of 



TTiR — rriQ — 2.5 log 



P 



D\2 / \ -2 / \ -2' 

n\ I r \ r — ri 



(6) 



where r is the particle's heliocentric distance, r\ — \ AU, vhq = —27.29 is the Sun's apparent 
R magnitude, Rmin = 20 km is adopted here, and the observation is presumed to occur at 
solar opposition. All of our calculations will also adopt the usual albedo of p = 0.04 so 
that our findings can be readily compared to past results obtained by others. However if 
an alternate albedo p is desired, simply revise all KBO radii reported here by a factor of 
(p/0.04)~^/^, and all masses by a factor of (p/0.04)~^/^. Finally, note that a power-law size 
distribution results in a cumulative luminosity function that varies as E(m) oc 10""*, where 
E(m) is the sky-plane number density of KBOs brighter than apparent magnitude m and 
the logarithmic slope is a — Q/5 (Irwin et al. 1995). 

Hubble Space Telescope observations reveal that the bright end of the Kuiper Belt's 
luminosity function has a steep logarithmic slope of et = 0.88 ± 0.1 for bodies having mag- 
nitudes niR -C 24, while the faint end {niR ^ 24) of the luminosity function has a shallow 
logarithmic slope oi a = 0.32 ± 0.15 (Bernstein et al. 2004); the steeper slope of the bright 
end of the luminosity function was also confirmed recently by Elliot et al. (2005). This lu- 
minosity function can be interpreted as evidence that the KBO size distribution is actually 
two power laws that break even at a magnitude rribreak — 24, which corresponds to a body 
of radius Rbreak — 65 km orbiting at a characteristic distance of r ~ 40 AU assuming it has 
an albedo of p = 0.04. However our apphcation will concentrate only on those KBOs that 
have known orbits, and 99% of those bodies have magnitudes rriR < rribreak- Consequently, 
this study will be sensitive only to the larger end of the KBO size spectrum, and such bod- 
ies will be characterized here via a single power-law size distribution having a — 0.88 and 
Q = 5q; = 4.4. 

Although ~ 5000 simulated particles in Fig. 3 manage to survive over the age of the 
Solar System, the Monte Carlo model assigns far too few of them with sizes that would 
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be detected by any telescopic survey of the Kuiper Belt. To boost the statistics of the 
detectable portion of the simulated population, each survivor is replicated 10^ times such 
that each particle's orbit elements a, e, i, Q, uj are preserved while its mean anomaly M is 
randomly distributed over < M < 2tt. It should be noted that this step effectively 
assumes that the particles' longitudes are uniformly distributed over 27r, which is not quite 
correct since Neptune's resonant perturbations tend to arrange the particles' longitudes in 
a non-uniform manner (Malhotra 1996; Chiang & Jordan 2002). Nonetheless, this is not a 
major concern since the observed KBOs were discovered along lines-of-sight that are roughly 
distributed uniformly in ecliptic longitude, which effectively washes-out Neptune's azimuthal 
arrangement of the Belt; see Appendix C for a more detailed examination and justification 
of this assumption. Lastly, this Monte Carlo model is then tested by verifying that the 
randomly generated population does indeed exhibit the expected luminosity function that 
varies as E(m) oc 10*5™/^. 

Further comparison of the Monte Carlo model of the Kuiper Belt to any observations 
must be done carefully. Note that the brighter KBOs tend to be discovered in shallow, 
wide-angle surveys that observe a large area Afi on the sky, while the fainter KBOs tend to 
be discovered in deeper surveys that observe smaller areas Af2. Consequently, the observed 
abundances of the various KBO subclasses (e.g., the Main Belt, the Scattered Disk, etc.) 
are proportional to all of these surveys' total area Jl(mR), which itself is some function of 
the limiting magnitude rriR. However Appendix B shows that this dependence upon Vt{mR) 
can be factored out by constructing ratios of the Belt's various subclasses. That Appendix 
also shows that the ratios of the observed abundance of any two dynamical classes of KBOs 
is approximately equal to the ratio of the intrinsic abundances of the much larger unseen 
populations^. Thus by plotting ratios of the simulated populations to the observed KBO 
populations, we can compare the model to the observations in a manner that is insensitive 
to survey details like their individual sky-coverage /STt. 

Figure 7 shows the apparent abundances of the 2:1 and the 3:2 populations relative 
to the Main Belt (MB) as a function of their i?-band magnitudes itlr- The upper light 
gray curve is the simulated ratio, which predicts an apparent 2:1 abundance that is about 
80% of the MB, while the dark gray curve is the observed ratio. Taking the ratio of these 
two curves reveals that their discrepancy at magnitudes uir > 21 (which refers to about 
90% of the observed sub-populations) is a factor of /2:i ~ 20 — the observed 2:1 resonance is 



^Of course, this method of analyzing the Belt's relative abundances will tells us little about those KBO 
populations that are either too rare, dim, or otherwise too difficult to recover in telescopic surveys. Nonethe- 
less, we still can use our method to place upper limits on the abundances of any hypothetical KBO populations 
that are unseen using the method described in Section 7. 



-14- 



markedly underabundant relative to the observed Main Belt population. There are Nmb{R > 
50 km) = 1.0 X 10^ Monte Carlo particles in the simulated Main Belt having radii R > 50 
km, and N2:i{R > 50 km) = 8.2 x 10"^ particles in the 2:1. If we let r2:i/MB represent the 
inferred ratio of 2:1 to Main Belt objects, then r2:i/MB — ^2:i//2:1^mb — 0.041, which is 
comparable (albeit lower by a factor of ~ 2) to the ratio that Trujillo et al. (2001) infer from 
telescopic surveys of the Kuiper Belt. 

The observed 3:2/MB ratio plotted in Figure 7 also shows that this resonant population 
is underabundant relative to model predictions by a factor of ~ 6 among bright objects 
with 21 < m/j < 23, and by a factor of ~ 60 at fainter magnitudes. Close inspection of 
the observations suggests that there indeed is a deficiency of fainter KBOs in the 3:2, and 
that this curve is not due to some overabundance of Main Belt KBOs having magnitudes of 
rriR ~ 23. 

It should be noted that the results given in Fig. 7 are not particularly sensitive to the 
detailed location of the Main Belt's outer edge. For instance, if we assumed the Belt's 
primordial edge where instead at a = 45 AU (e.g., Trujillo & Brown 2001; see also Section 
7), this reduces the 2:1 and MB populations both by about 40% while leaving the 3:2 
population unchanged. Consequently, the 2:1/MB ratios of Fig. 7 are largely unchanged for 
both the simulated and observed populations, while the 3:2/MB ratios increase by a factor of 
(1 — 0.4)^^ ~ 1.7. However the discrepancy between the simulated and observed populations 
is still the same factors of ~ 6-60. 

Although there are several possible interpretations of the discrepancies seen Fig. 7, the 
most plausible explanation is that other unmodeled processes are responsible for (i) reducing 
the trapping efficiencies of the 2:1 and 3:2 resonances by factors of ~ 6-60, or (m.) causing 
trapped particles to diffuse out of the resonances and into nearby regions of phase space 
that are quite unstable (cf. Fig. 1 of Duncan et al. 1995), resulting in their ejection from the 
Kuiper Belt. Such unmodeled processes include the collisions and gravitational scatterings 
that occurred with ever greater vigor during earlier times when the Belt was more crowded. 
The scattering of these planetesimals by Neptune was of course responsible for driving that 
planet's migration, so the occasional scattering of a large and/or close planetesimal will 
cause that planet's orbit and hence its resonances to shudder some. Likewise, scattering 
events among the KBOs themselves would also cause their semimajor axes to diffuse some, 
as would collisions. This means that scatterings and collisions will have driven a random 
walk in the resonant particles' semimajor axes, as well as a random walk in the location 
of the resonances themselves. It is possible then that these unmodeled effects can drive 
particles out of resonances and reduce the resonant population by the large factors indicated 
by Fig. 7, a scenario that is also explored in simulations by Zhou et al. (2002) and Tiscareno 



- 15 - 



& Malhotra (2004). 

The magnitude dependence of the observed 3:2/MB ratio shown in Fig. 7 is also quite 
curious. The fact that the observed ratio varies with apparent magnitude mn, while the 
simulated ratio remains constant at magnitudes fainter than niR — 21, suggests that the 
N{R) (X power-law that was universally applied throughout the entire Belt is overly 
simplistic. One way for the model to achieve better agreement with the observations is to 
assume that larger, brighter bodies are more abundant in the 3:2, and that smaller, fainter 
bodies are less abundant there than they arc in the MB, which requires a shallower size 
distribution. The dashed curve in Fig. 7 illustrates this possibility, which shows the simulated 
3:2/MB ratio assuming that the 3:2 bodies have a shallow Q = 2.7 size distribution with 
Rmin — 4.3 km (note that reducing Rmin has the effect of reducing the total number of visible 
objects) while the MB bodies have the usual distribution with Q — 4.4 and Rmin — 20 km. 
The shallow size distribution that is inferred here for the 3:2 population is also consistent with 
the logarithmic slope of a ~ 0.56 that Elliot ct al. (2005) recently reported for the luminosity 
function of their 'resonant' population that is dominated by 3:2 KBOs; the size distribution 
inferred from that work is Q = 5q; = 2.8. KBO sizes can also vary with inclination [Levison 
& Stern (2001), but see also footnote 1 of Gomes (2003)]. In particular, Bernstein et al. 
(2004) report that the bright end of the luminosity function for high-inclination [i > 5°) 
KBOs have a logarithmic slope of a — 0.66 and a size distribution Q — 5a — 3.3 that is 
much shallower than the low-i KBOs having a = 1.36 and a Q = 5a — 6.8. 

In our Monte Carlo model there arc only N^.2{R > 50 km) = 2100 bodies larger than 
R = 50 km, so their numerical abundance relative to the Main Belt is r^a/MB = N^;2/Nmb = 
0.021, which again is comparable (but again lower by a factor of ~ 2) to the ratio reported 
in Trujillo et al. (2001). The largest Monte Carlo body at the 3:2 has a radius R ~ 1000 km, 
which is comparable to the size of the largest multi-opposition KBO there^ with R ~ 1100 
km, assuming p = 0.04. 

This range of power-law indices that is inferred for the Kuiper Belt, 2.7 < Q < 4.4, 
is comparable to the values of Q that are observed at various sites throughout the asteroid 
belt. The Near Earth Objects have a fairly shallow size distribution with Q — 1.95 (Stuart 
& Binzel 2004), while the asteroid families exhibit steeper size distributions. For instance. 
Fig. 1. of Tanga et al. (1999) shows size distributions for several prominent asteroid families 
having values of 2 < Q ^ 6. Note also that nonfamily asteroids have Q ~ 3.0 (Ivezic 
et al. 2001), which is slightly steeper than the canonical Q = 2.5 value that results from 
a coUisional cascade (Dohnanyi 1969). Since the various asteroid subclasses exhibit such a 



excepting Pluto of course which has an anomalously high albedo of p ~ 0.5. 
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wide variation in their size distributions over a relatively narrow range of semimajor axes of 
Aa ~ 4 AU, perhaps it should be of no surprise that the spatially much wider Kuiper Belt 
might also exhibit some variety in Q. 

7. The outer edge of the Solar System 

Inspection of Figure 3 shows a prominent absence of observed KBOs having modest 
eccentricities of e ~ 0.1 near and beyond Neptune's 2:1 resonance. The prevailing interpre- 
tation of this observed feature is that there is a boundary near 45 < a < 50 AU that marks 
the outer edge of the Solar System's primordial Kuiper Belt (Allen et al. 2001; Trujillo & 
Brown 2001). The dots in Fig. 8 shows the Belt's surface density o"(r) that is inferred from 
these observations, which peaks at r ~ 45 AU. We have nonetheless allowed our simulated 
Kuiper Belt to extend out to a = 80 AU in order to use the dearth of observed distant KBOs 
to place quantitative upper limits on the abundance of hypothetical KBOs that might live 
beyond 50 AU. 

The Nbody/Monte Carlo of the previous Section can be used to predict how many KBOs 
should have been observed in the a > 50 AU zone (which we identify here as the Outer Belt, or 
OB) assuming (i.) that the primordial Kuiper Belt extends smoothly out to a = 80 AU, and 
(ii.) that all KBOs everywhere have the same size distribution with the usual parameters 
Q — 4.4 and Rmin — 20 km. This simulation's ratio of Outer Belt/Main Belt objects, 
t^ob/mb{''^ir)- is plotted versus magnitude uir in Fig. 9. This is the ratio of Nq'^'Ittir) — 
the number of bodies in the simulated Outer Belt (whose members have semimajor axes 
50 < a < 80 AU and eccentricities e < 0.2), to N'^^{mfi) = the number of bodies in the 
simulated Main Belt (where 40.1 < a < 47.2 AU) in the magnitude interval m^i Am where 
Am = 0.5. According to the Figure, the expected OB/MB ratio is tob-.mb = 0.4. At present 
there are A^mb — 264 KBOs in the MB that have been observed for 2 or more oppositions, 
and the dimmest member of this group of KBOs has an apparent magnitude = 24.5. The 
Nbody/Monte Carlo model thus predicts that there should also be Nq^ — vob/mb^mb — 
100 objects brighter than orbiting in the OB beyond a = 50 AU. This prediction is in 
marked contrast with the observations which show that there are no known multi-opposition 
objects orbiting in the OB with magnitudes brighter than m^, i.e., Nq^ < 1. 

This prediction that the Outer Belt would have an observed abundance that is 40% of 
the Main Belt differs considerably from that of Gladman et al. (1998) who estimated that 
the observed Outer Belt population should only be ~ 6% of the total observed population. 
However this much lower estimate was obtained by assuming that the current Belt's surface 
number density resembles its primordial a{r) which likely varied as or so. This very 
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common assumption causes the inner part of the Belt to be more concentrated than its 
outer part. However a real Kuiper Belt has been dynamically eroded over the eons by the 
giant planets' gravitational perturbations. This dynamical erosion is illustrated in Fig. 8, 
which shows the simulated Belt's primordial surface density (gray curve) and its final surface 
density (dashed curve); similar erosion is also seen in the long-term integrations of Duncan 
et al. (1995). Figure 8 shows that a{r) for an eroded Belt is a sharply increasing function 
of r for r < 50 AU, which imphes that the inner observable portion of the Belt is very 
underdense relative to the more distant r ~ 50 AU zone. This dynamical erosion accounts 
for the discrepancies between our Outer Belt predictions and that by Gladman et al. (1998). 

Recall that Figs. 3 and 7 show that the 2:1 and 3:2 populations are very depleted relative 
to model predictions, and that the zone beyond Neptune's 2:1 resonance is either empty or 
inhabited by bodies too small and faint to be seen. To account for these depletions, the solid 
curve in Fig. 8 also shows a revised surface density curve that is obtained from the simulated 
Belt that is truncated at a = 45 AU (about 3 AU inwards of Neptune's 2:1 resonance), and 
with the negligible contribution from the 3:2 populations also being ignored. This result is a 
curve that agrees quite well with the Belt's observed surface density variations. Despite this 
good agreement in the radial distributions of the simulated and observed Kuiper Belts, Fig. 
3 shows that this apparent edge at a = 45 AU is still rather fuzzy since there are four multi- 
opposition KBOs of low eccentricity (e < 0.05) orbiting in the Main Belt at 45 < a < 48 
AU. Close inspection of Fig. 3 shows that a hard edge at a = 45 AU also could not account 
for the KBOs having e ~ 0.1 in the 45 < a < 48 AU, unless the advancing 2:1 resonance 
also dragged some bodies out of the a < 45 AU zone and deposited them here, reminiscent 
of the scenario suggested by Levison & Morbidelli (2003). 

The remainder of this Section places upper limits on the size and abundance of any 
unseen KBOs that might lurk beyond a = 50 AU. Of course there are multiple interpretations 
of the dearth of observed multiple-opposition bodies orbiting beyond a = 50 with modest 
eccentricities of e ~ 0.1, i.e., that Nq^ < 1. One interpretation of this upper limit is 
that assumption (i) is incorrect — that the primordial Kuiper Belt's density did not extend 
smoothly beyond Neptune's 2:1, but that it instead was reduced by a factor / (relative to 
the smooth model's density) in the a > 50 AU zone. In this case, the OB/MB ratio becomes 
roB/MB = 0.4// = N^^^/N^% < l/N^%, which implies that the primordial density of the 
OB was smaller than the MB by a factor / > 100. 

Alternatively, assumption (ii.) could be incorrect, namely, that the KBO size distribu- 
tion was not uniform everywhere. For instance, the absence of any multi-opposition bodies 
in the OB having magnitudes brighter than = 24.5 could simply mean that bodies be- 
yond r ~ 50 AU are dimmer than and thus have radii smaller than R ~ 80(p/0.04)~^/^ 
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km (see Eqn. 6). Note that Trujillo et al. (2001) obtained a similar limit, but that they came 
to regard this scenario as unlikely. 

It is also possible that the Outer Belt's size distribution is steeper, i.e., has a larger 
Q, than the Main Belt's size distribution. An increase in Q decreases the abundance of 
bright bodies, as is illustrated by the curve in Fig. 9 which gives the the tob/mb ratio for 
an OB having a, Q = 6.0 size distribution while bodies in the MB have the usual Q = 4.4 
distribution. This particular Q is also the minimum value that is consistent with the observed 
upper limit of Tob/mb < 1/264; Outer Belts with with a smaller Q would contain at least 
1 KBO brighter than m*^ in the a > 50 AU zone for every 264 KBOs detected in the Main 
Belt, while an OB having a larger Q would be undetected. This particular model is near the 
threshold of detection, and its largest member has a radius oi R — 250 km. Trujillo et al. 
(2001) also considered this scenario, but they concluded that the absence of distant KBOs 
requires a steeper Q > 9 size distribution. The origin of this discrepancy is unclear. 

It is also interesting to note that the low-inchnation KBOs have a logarithmic slope of 
a — 1.36 along the bright end of their luminosity function (Bernstein et al. 2004), which 
implies a steep size distribution of Q = 5q; = 6.8. Such bodies, if they inhabit the Outer 
Belt beyond a = 50 AU with the same abundances as adopted by our model, could conceiv- 
ably have avoided detection to date due to their steep size distribution. In other words, a 
Main Belt whose low-i population extends beyond a = 50 AU while its high-i population 
terminates at a = 50 AU could be quite consistent with their non-detection. 

Prom these considerations it may be concluded that the observed absence of multi- 
opposition KBOs in the a > 50 AU zone having modest eccentricities e ~ 0.1 implies: {i.) 
that this part of the primordial Kuiper Belt was underdense by a factor of / > 100 relative 
to the a < 50 AU zone, or that (m.) these distant KBOs have radii i? < 80 km, or that (m.) 
their size distribution has a power-law index Q > 6.0, or perhaps (iv.) some combination of 
the above effects. 



8. The origin of Centaurs 

It is generally accepted that Centaurs are those bodies that have diffused inwards from 
the Kuiper Belt into orbits that cross the giant planets (e.g., Duncan et al. 1988). Presently 
there are 27 known Centaurs observed for more than one opposition; these are the red dots 
in Fig. 10 having a < 30 AU. Since planet-crossers are quickly ejected or accreted. Centaurs 
have short dynamical hfetimes of only ~ 10^ years (Levison & Duncan 1997; Tiscareno & 
Malhotra 2003). Consequently, the density of these 'escapees from the Kuiper Belt' (e.g.. 
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Stern & Campins 1996) is very tenuous inside of a = 30 AU (see Fig. 8). Indeed, only 
Nq — 7 Centaurs are detected during the final AT = 2 billion years of our simulation that 
was only sparsely time-sampled once every At = 10® years, so the instantaneous number of 
Centaurs is Nq = N^lAt / AT) = 0.35 at the end of the simulation. There are also Nmb = 
565 bodies in the Main Belt, so the Centaur/Main Belt ratio is provisionally estimated at 
rc/MB = 6.2 X 10-"^. 

The open circles on Fig. 10 show the orbital elements of these seven Centaurs at time 
t = 10® years, which is at a time when planet migration has only recently ceased. Thus 
the open dots indicate the locations where Neptune has parked these proto-Ccntaurs in the 
Kuiper Belt. Note that all seven Centaurs originate from sites in/near Neptune's mean- 
motion resonances, namely, the 3:2, 5:3, 13:7, 2:1, and the 5:2. Their subsequent motions at 
times t > 10® years are shown as black dots (again, poorly time-sampled), which show that 
the eccentricities of nearly all proto-Centaurs initially wander up-and-down with constant 
a until they have a close encounter with Neptune, scatter off that planet, make a brief 
apparition in the a < 30 AU Centaur zone, and then are quickly removed from the system. 

These seven bodies have initial semimajor axes of 28 < a < 48 AU at time t — 0, so 
Centaurs can also be regarded as samples that have been drawn from a wide swath of the 
outer solar nebula. Fig. 10 also shows that the simulated Centaurs are concentrated just 
inside of Neptune's orbit; their mean heliocentric distance is r = 26 ± 3 AU, and their mean 
inclination is ?' = 16 ± 10°. Note also that three of the seven Centaurs emerged from the 3:2 
and 2:1 resonances, which Section 6 showed to be heavily depleted relative to the model's 
predictions. Consequently, the Centaur/Main Belt ratio reported above should instead be 
interpreted as an upper limit, e.g., tc/mb < 6.2 x 10~^. It is shown later in Section 11 that 
our model predicts that there are Nmb — 1-3 x 10^ Main Belt KBOs having radii it! > 50 km, 
so this model also predicts that there are Nq — tc/mbNmb < 80 similarly-sized Centaurs. 

Although the Centaur upper limit reported here is comparable to the population that 
Sheppard et al. (2000) infer from the Centaur luminosity function, there is still a promi- 
nent disconnect in the heliocentric distances of the simulated and observed populations; our 
simulated Centaurs all reside at r > 22 AU, while the three Centaurs that Sheppard et al. 
(2000) used to construct the Centaur luminosity were detected at heliocentric distances of 
r < 19 AU. One possible interpretation of this excess of Centaurs at r < 20 AU is that 
Centaurs may be breaking up and spawning new Centaurs (e.g., Pittich & Rickman 1994) as 
they wander among the giant planets. Finally, we note that deep, wide-angle surveys of the 
Kuiper Belt, such as the Legacy Survey that is currently being implemented at the Canada 
France Hawaii Telescope, may soon reveal the existence of the more distant Centaurs that 
is anticipated by this model to reside at distances of 23 < r < 29 AU. 
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9. Neptune's Trojans 

Figure 3 also shows that Nt = 5 particles managed to survive the length of the simula- 
tion at Neptune's 1:1 resonance. These simulated particles are of course Neptune's Trojans, 
of which two are presently known: 2001 QR322 (Chiang et al. 2003) and 2004 UPio (Sheppard 
et al. 2005). For this simulation the Trojan/MB ratio is Tt-.mb = Nt/Nmb = 8.8 x 10~^, 
where Nmb — 565 is the number of survivors that persist in the Main Belt. The spatial 
coordinates of the two observed and five simulated Trojans are shown in Figure 11, which 
indicates that these particles can roam about with longitudes of A0 ^ ±30° from Neptune's 
triangular Lagrange points with semimajor axes of Aa ~ ±0.32 AU from Neptune's. The 
extent of these Trojan sites are similar to that seen in integrations by Holman & Wisdom 
(1993) and Nesvorny & Dones (2002). Note that no special effort was made to start any of 
the simulated particles at Neptune's Lagrange points. Rather, all particles were distributed 
randomly about a disk according to a smooth surface density law, with the inner edge of 
the disk being well inside of Neptune's initial tadpole region. In our simulation the five 
survivors had initial semimajor axes of Aa ~ ±0.28 AU from Neptune's initial a, and there 
were a total of 68 particles initially in Neptune's Trojan source region (i.e., |Aa| < 0.28 
AU and |A0| < 30°), so the surviving Trojan fraction is about 7%. This survival fraction 
is comparable to that obtained by Kortenkamp et al. (2004) in a similar simulation. That 
work also showed that as planets migrate, several secondary resonances sweep across the 1:1, 
which result in a heavy loss of Neptune's Trojans during the migration epoch. 

Neptune's Trojans are of interest since they might place constraints on some models of 
the early evolution of the outer Solar System. For example, Thommes et al. (1999, 2002) 
postulate that Neptune originally formed in the vicinity of Jupiter and Saturn and was tossed 
outwards after scattering off the larger planets. But the existence of 2001 QR322 and 2004 
UPio might cast doubt on this scenario since Trojans seem unlikely to persevere at Neptune's 
Langrange during such a scattering event. However it has since been shown that a recently- 
scattered Neptune can still acquire its Trojans later as that planet's orbit is circularized by 
a dense Kuiper Belt (Levison 2005, personal communication). It is also conceivable that 
Neptune may have captured its Trojan from a heliocentric orbit after Neptune's orbit has 
settled down. Although Kortenkamp et al. (2004) shows that the direct capture of Trojans 
from heliocentric space is rare and results in only transient Trojans, Chiang & Lithwick 
(2005) show that mutual collisions may have inserted small bodies into stable orbits at 
Neptune's Lagrange points after its orbit has circularized. 
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10. The Extended Scattered Disk 

Figure 12 shows the orbits of those scattered particles that have been tossed into very 
wide orbits about the Sun. Most of the simulated scattered particles have perihelia between 
30 < q' < 40, as do most of the observed scattered KBOs. However there are two exceptions 
to this rule, namely, 2000 CR105 and 2003 VB12 (also known as Sedna) which have respective 
periheha of g = 44.14 ± 0.02 (Gladman et al. 2002) and 5 = 76 ± 4 AU (Brown et al. 2004). 
Gladman et al. (2002) classify those Scattered KBOs having periheha higher than g ~ 40 
AU as members of a so-called extended Scattered Disk. Sedna's large radius of ~ 1000 
km makes this object a particular curiosity since its discovery circumstance suggests that 
there may be a few hundred other unseen Senda-sized objects (Brown et al. 2004). Since 
Sedna has a mass of ~ 10~^ M^, the implied mass that might be hidden in the extended 
Scattered Disk is a few tenths of an Earth-mass. Thus Sedna by itself may represent an 
enormous reservoir of unseen mass that is comparable to the 'conventional' Kuiper Belt (see 
Section 11). 

The extended Scattered Disk is also of dynamical interest since, as Gladman et al. 
(2002) note, most dynamical models of the Kuiper Belt (including this one) generally produce 
Scattered Objects in lower periheha orbits having g < 40 AU. Gladman et al. (2002) review 
a number of scenarios that might explain how a KBO might get promoted from a nearly 
circular orbit into a wide, eccentric orbit having ? > 40 AU; these include: (i.) chaotic 
diffusion of scattered bodies, (ii.) gravitational scattering by long-gone massive protoplanets, 
{Hi.) scattering by an undiscovered distant planet, and (iv.) scattering by a single star that 
passes to within ~ 100 AU of the Sun. However all of these scenarios are problematic. For 
instance, billion-year integrations of ~ 10^ particles in chaotic Neptune-scattered orbits fail 
to diffuse into orbits having perihelia as high as that of 2000 CR105 (Gladman et al. 2002). 
Morbidelli et al. (2002) also cast doubt on scenarios {ii-iii.) by showing that ~ 20% of any 
distant population of protoplanets would have persisted over the age of the Solar System, 
and that some fraction of these large bodies should already have been discovered by one of 
the various wide-angle Kuiper Belt surveys. Scenario (iv.) is also in doubt since simulations 
of a close encounter with a single star generally produce disturbances in an outer Kuiper 
Belt that is quite unlike that seen in the observed Belt (e.g., Ida et al. 2000). However 
Fernandez & Brunini (2000) have shown that repeated encounters by more distant stars can 
produce Sedna-like orbits. This may have occurred early while the Sun was still a member 
of the open cluster from which it presumably formed. In this scenario, the giant planets 
scatter small bodies into wide orbits of a ~ 100 to 1000 AU, which the nearby cluster stars 
then perturb into Sedna-like orbits having higher perihelia. Scattering by a passing star was 
recently re-examined by Morbidelli & Levison (2004), and their simulations also support 
this scenario. 
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Although our simulations did not produce any Sedna-like objects in orbits that are well- 
decoupled from the giant planets, wc did find a single scattered object in a 2000 CRio5-like 
orbit in the Extended Scattered Disk with a semimajor axis of a = 92 AU and a perihelion of 
q = 43.5 AU (the large black dot in Figure 12). The orbital history of this scattered particle 
is shown in Fig. 13, which shows that as the particle inhabited Neptune's 16:3 resonance 
during times 2 < t < 3 Gyrs, some process raised this Scattered particle's perihelion up and 
into the Extended Scattered Disk on a biUion-year timescale. This kind of behavior was 
first reported in Levison & Duncan (1997) and Duncan & Levison (1997), whose simulations 
also show that some scattered particles can achieve high perihelia orbits while in or near 
mean-motion resonances. 

However we have not identified any particular resonance as being responsible for raising 
the perihelia of our one CRios-candidate shown in Fig. 13. For instance, a Kozai resonance 
is not implicated since the argument of perihelion cu does not librate. The possibility of other 
Pluto-like 'super-resonances' (e.g., Malhotra & Williams 1997) was also examined; this is the 
libration of a resonance angle of the form = j {uj — ujn) — k{Q — Q n) , where Q is the particle's 
longitude of the ascending node and the subscript refers to Neptune's orbit elements. 
Angles having < \k\ < 10 were examined, and although the angle = 5(f2 — Q^) angle 
did in fact librate for about 1 Gyrs, that occurred well after the time when the particle's q 
was raised. A resonance involving interactions with multiple planets is also unlikely since 
the particle's Tisserand parameter T (which is simply its Jacobi integral sans the interaction 
energy due to Neptune) was well-preserved during these times. Although the particular 
mechanism that drove this particle into the Extended Scattered Disk in not understood, this 
particle docs demonstrate that it is indeed possible for Scattered particles to diffuse into the 
Extended Scattered Disk via planetary perturbations alone, with other external agents (like 
a stellar encounters) being absent. This transport from the Scattered Disk to the Extended 
Scattered Disk via mean-motion resonances is also evident in the simulations of Levison & 
Duncan (1997) and Duncan & Levison (1997). However this transport has an extremely low 
fiux since only one of the ~ 5 x 10^ particles initially in the Scattered Disk did manage to 
enter the Extended Scattered Disk and persist over the age of the Solar System. External 
perturbations from passing stars (Fernandez & Brunini 2000; Morbidelli & Levison 2004) 
may indeed be more effective at producing members of the Extended Scattered Disk. 

Note also the simulated particles represented by crosses in Fig. 12. Even though their 
perihelia of 42 < g < 54 AU might suggest that they also inhabit the domain of the Extended 
Scattered Disk, they are in fact resonant particles that were trapped at the 3:1, 7:2, and 4:1 
resonances during the migration epoch. Most of these particles have libration amplitudes 
less than |A0jfc| < 90°. These particles also had initial semimajor axes of a > 47 AU, 
which is noteworthy since, if any resonant KBOs are ever discovered in these orbits, they 
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could be interpreted as evidence that the outer edge of the Solar System hes beyond a > 47 
AU. However that interpretation would still be ambiguous, since Neptune-scattered evaders, 
which originated from smaller semimajor axes, can also settle into these same resonances — 
see Fig. 18 of Gomes (2003) for an example. 

10.1. The Scattered Disk 

Figure 14 shows the apparent abundance of so-called Scattered Disk objects relative to 

the Main Belt, as predicted by the Nbody/Monte Carlo model. There are Nsd = 1-9 x 10^ 
Monte Carlo bodies in orbits having 50 < a < 150 AU and perihelia 28 < g < 40 AU, while 
Nmb = 1-0 X 10^ Monte Carlo particles survive in the Main Belt, so the model predicts 
an intrinsic SD/MB ratio of Tsd/mb = Nsd/Nmb = 0.19. The apparent ratio of these two 
populations is ~ 0.1, also shown in Fig. 14. Note that the intrinsic SD/MB ratio inferred 
here is about one-fourth that reported in Trujillo et al. (2001). 

11. Calibration 

The ecliptic luminosity function of Bernstein et al. (2004) is shown in Fig. 15, and 
its bright end varies as T,(mji) = 10"('"^™") deg~^ where a = 0.88 and tuq = 23.1, for 
magnitudes ijir < rrihreak = 24. This luminosity function gives the number density of KBOs 
near the ecliptic that are brighter than magnitude ttir. Since this curve scales with the total 
number of KBOs, it can be used to calibrate the simulation to determine the total number 
of objects in the Kuiper Belt. 

Sections 6-7 show that the observed 3:2 and 2:1 populations are severely depleted rela- 
tive to model predictions, and that bodies in an Outer Belt beyond the 2:1 are either absent 
or too faint to be seen. To account for these depletions, the truncated Kuiper Belt similar 
to that of Section 7 is adopted; this Belt is formed by discarding any bodies orbiting beyond 
the 2:1, as well as all bodies orbiting within Aa = 0.6 AU of Neptune's 3:2 and 2:1 reso- 
nances. There are Nj^i^ody = 587 Nbody particles in this truncated Kuiper Belt, and they 
are replicated A^reps = 10^ times with sizes and magnitudes assigned to them according to 
the Monte Carlo method of Section 6, with Q = 4.4 and Rmin = 20 km. 

The simulation's median inclination is low (e.g.. Section 5), only isim = 2.7°, which 
is much lower than the median inclination iobs = 15.6° that is inferred from the debiased 
KBO inclination distribution reported by Brown (2001). Due to these low inclinations, 
the simulation's ecliptic luminosity function would thus be artificially overdense by a factor 
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fi — iohs/isim — 5.8, SO it is revised downwards by this factor to compensate. The simulation's 
E(mR) is then multiphed by a factor /^r = 3.5 to fit it to the bright end of the observed 
luminosity function; this accounts for the different populations in the simulated and observed 
Kuiper Belts, and results in the curve shown in Fig. 15. The size distribution adopted here 
is valid down to a radius of about Rbreak = 65 km (see Section 6), so the inferred number 

of KBOs larger than Rbreak is N^reak = NmodyNrepsf N {Rbreak / Rminy^ - LI X 10^. To 

estimate the total number of KBOs larger than the fiducial radius of R50 — 50 km, note 
that the faint end of the observed luminosity function has a logarithmic slope of a faint — 
0.32 (Bernstein et al. 2004), which implies a power-law index of Qfaint — ^OLfaint — 1.6 
for bodies having radii R < Rbreak- The total number of bodies larger than R^q is thus 

A^50 = Nbreak{R50/Rbreak)-'^''''-' ^ 1.7 X 10^ 

The total mass of these bodies is obtained from their cumulative size distribution, which 
for the large bodies with it! > Rbreak can be written N{R) — Nbreak{R/ Rbreak)~'^ ■ The 
differential size distribution is then dN{R) = \dN/dR\dR, and if M(R) = mass of a body 
having a radius R, the total mass of bodies having radii in the interval Rbreak < R < Rmax 
is 



M{R> Rbreak) 
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break 
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NbreakMbreak (7) 



where Mbreak = M{Rbreak) — 1-2 X 10^^(p/l gm/cm^)(p/0.04)~^/^ gm, which is the mass of 
a body of radius Rbreak = 65 km assuming it has a density p and albedo p. The total mass 
of KBOs larger than Rbreak with semimajor axes inside of Neptune's 2:1 is thus M{R > 
Rbreak) — 0.07(p/l gm/cm^)(p/0.04)~^/^ for a Q = 4.4 size distribution that extends 
to radii as large as Rmax — 1000 km. To get the total mass of bodies at the fiducial size 
R — i?5o, add to the above the mass of bodies in the size interval -R50 < R < Rbreak, which 
is roughly AM ~ (A^5o — N break) Mbreak- The total mass of bodies larger than i?5o = 50 km 
is then Mtotai = M{R > Rbreak) + AM ~ 0.08(p/l gm/cm3)(p/0.04)-3/2 M®. Note that the 
0.08 prefactor is a consequence of adopting the oft-employed Halley albedo oi p = 0.04. 
However recent observations indicate KBOs have an average albedo of p ~ 0.1 (Altenhoff 
et al. 2004; Grundy et al. 2005), which in turn lowers the Kuiper Belt mass to Mtotai 0.02 
Mq assuming they have a unit density. 

This population estimate is comparable to, but a bit higher than, previous estimates 
that rely on far simpler models of the Kuiper Belt. For instance, Trujillo et al. (2001) 
report a Main Belt population of 3.8 x 10^ objects of mass 0.03 among bodies having 
radii it! > 50 km. They also estimate the Belt's total population to be 1.9x the Main Belt 
population, so a total population of A^5o ~ 7.2 x 10^ bodies larger than it!5o having mass 
Mtotai — 0.06 M0 is inferred. A similar estimate is also inferred from the HST survey by 
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Bernstein et al. (2004); according to their Fig. 8, the sky-plane number density of KBOs 
larger than R^q is E(i? > /?5o) ~ 13 deg~^. Since the Kuiper Belt subtends a total solid 
angle of AQfotai — 8100 deg^ (Brown 2001), the total number of KBOs larger than R^q is 
A^50 = E(i? > R5o)A9,totai — 1-1 X 10^ having a total mass of Mtotai — 0.05 M®. 

Sections 6-10 show that the simulated Belt's various dynamical classes have abundances 

of r2:i/MB = 0.041, r3;2/Afij = 0.021, rc/MB < 6.2 x 10""^, tt/mb < 8.8 x lO'^, and tsd/mb = 
0.19 relative to the Main Belt, so the Main Belt fraction is J'mb = 1 — '^i'^^i/MB — 0.74 
and thus there are are Nmb{R > 50 km) = /mbN^q ~ 1.3 x 10^ Main Belt KBOs having 
radii i? > 50 km. The numerical abundance of the i^^ dynamical class is Ni{R > 50 km) = 
Ti/MBfMBN^o, and its mass is Mi{R > 50 km) = Ti/MBfMBMtotai where Mtotai = 0.08 M® 
assuming p = I gm/cm^ and p — 0.04; these abundances and masses are listed in Table 1. 
The exception is the 3:2 mass estimate which adopts the Q — 2.7 power-law size distribution 
described in Section 6; if this subgroup really does have such a flat size distribution, then 
Eqn. (7) must be used to calculate its mass^. 

Note also that the preceding ratios assume that the Main Belt terminates just inwards 
of the 2:1 resonance at a = 47.2 AU. If, however, one wishes to adopt an outer edge at a = 45 
AU, then Section 6 shows that this reduces the MB population by 40%, so that the ratios 
r quoted above should then be raised by a factor of 1.7. The exception to this rule are the 
bodies at the 2:1 resonance — their abundance relative to the MB is unchanged. However the 
total number of KBOs reported here is still insensitive to the detailed location of the Main 
Belt's outer edge, since that number is obtained by fitting the simulated KBOs' luminosity 
function to the observed S(mij), which is quite insensitive to the detailed location of the 
Main Belt's outer edge. 

It should also be noted that this study employed an initial cr{a) oc disk surface den- 
sity, but our findings are readily adapted for an alternate surface density law. For instance, 
if the canonical a{a) oc a^^'^ law were instead desired, then this shallower power law would 
result fewer objects trapped in the 3:2 resonance relative to the Main Belt population. Since 
the 3:2 objects are drawn from the 03:2 ~ 32 AU part of the disk, while the Main Belt objects 
formed at qmb ~ 44 AU, this revised surface density law would reduce the 3:2/MB ratio 
reported here by a factor (a3;2/aMB)°'^ ~ 0.85, which is a 15% change in relative abundance. 
Of course, the 2:1/MB ratio would remain unchanged since their source populations are the 
same. 

Also, we conservatively interpret the abundance of Neptune's Trojans reported Table 1 
as an upper limit on their real abundance. It was argued in Section 6 that other unmod- 



^with the quantities Rbreak, Nfyreak, and Mbreak replaced by R5Q, r3:2/MBNMB, and M{R = 50 km). 
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eled processes, possibly the scattering of planetesimals by Neptune or amongst themselves, 
reduced the trapping efficiency of the 3:2 and 2:1 resonances by factors of ~ 10. Thus it is 
possible that the same unmodeled phenomena might also have destabilized orbits at Nep- 
tune's 1:1 resonance, so the actual number of Trojan survivors may be smaller than that 
reported in Section 11. 

Finally, upper limits on the abundance of KBOs inhabiting a hypothetical Outer Belt 
are reported for the 50 < a < 80 AU zone assuming these bodies have the shallowest possible 
size distribution, namely Q = 6.0 down to Rmin = 20 km (see Section 7). In this case, there 
are at most Nqb = 1-3 x 10^ bodies in the OB having radii of 50 < /? < 250 km and a total 
mass of Mob — 0.008 assuming a density of p = 1 gm/cm'^ and an albedo p = 0.04. 

12. Effects not modeled 

It should be noted that the model used here only accounts for the Belt's dynamical 
erosion that is a consequence of Neptune's gravitational perturbations — it does not account 
for the collisional erosion of the Kuipcr Belt that is often invoked to account for the Belt's 
depleted appearance (e.g.. Stern 1996; Kenyon & Luu 1999). In particular, models of KBO 
accretion, as well as the self-consistent Nbody simulations of Neptune's migration, all suggest 
that the Kuiper Belt's primordial mass was of order ~ 30 (Stern 1996; Kenyon & Luu 
1999; Hahn & Malhotra 1999; Gomes et al. 2004), which is at least ~ 400 times more than the 
current mass. However the model used here, which only accounts for the dynamical erosion, 
results in a depletion by a factor of about 3 in the 30 < a < 48 AU zone of Fig. 3. This 
suggests that collisional erosion, which is not modeled here, may have been responsible for 
reducing the Belt's mass by an additional factor^ of ~ 100. Nonetheless, the abundance and 
mass estimates obtained here should still be reliable provided the Belt's collisional erosion 
was relatively uniform across the observable 35 < a < 50 AU zone. If, however, collisional 
erosion was more vigorous in some parts of the Belt, and less so in other parts, then the 
estimates obtained above will only be accurate in an order-of-magnitude sense. 

A comparable problem also occurs with the model's inclinations. Section 11 shows that 
the simulated Kuiper Belt is too thin by a factor of fi ~ 6. This is compensated for by 
reducing the simulation's luminosity function E(mi^) by the factor fi, which is equivalent to 
increasing each particles' inclination by this factor. Again, this crude treatment should still 



^Only a 'conventional' Kuiper Belt model, like the one explored here, need invoke additional erosion to 

reduce the Kuiper Belt mass by another factor of ^ 100. This is distinct from the push-out model which 
need not rely on any collisional depletion of the Kuiper Belt (Levison & Morbidelli 2003). 
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yield a reliable estimate of the KBO population provided the factor /j is uniform everywhere 
and independent of semimajor axis a. If, however, /j is not independent of a, then this will 
result in errors in the relative abundances of the Belt's various subpopulations reported in 
Table 1. 

We also note that the relative abundances of the Belt's various subpopulations are de- 
termined by a model that invokes a smooth outward migration by Neptune by Aajv — 9 AU, 
with the results reported in Fig. 6. That Figure shows that the smooth migration scenario 
predicts a combined 3:2 + 2:1 population that is comparable to the Main Belt population. 
This is because smooth migration is very efficient at trapping particles at Neptune's reso- 
nances, and this results in densely populated resonances. However a detailed comparison 
of the model predictions to the observed abundances indicates that the resonant KBO pop- 
ulation is really only about 5% of the Main Belt population (see Section 6 and Table 1). 
The seemingly low abundance of resonant KBOs is hkely due to unmodeled effects that may 
have occurred during the migration epoch, possibly due to the mutual scattering that might 
occur among bodies trapped at resonance, or perhaps due to the gravitational scattering of 
large planetesimals by Neptune (e.g., Zhou et al. 2002; Tiscareno & Malhotra 2004). 

There is also evidence indicating that a wide swath of the early Kuiper Belt was stirred 

up prior to the onset of Neptune's migration. Recall that simulations of Neptune's outward 
migration into a dynamically cold Kuiper Belt is unable to account for the eccentricities of 
e ~ 0.1 observed among Main Belt KBOs (Section 2). This suggests that the Belt was stirred 
up, cither prior to or after the onset of Neptune's migration. However Section 3 shows that 
this stirring event likely occurred prior to migration: migration into a stirred-up Kuiper Belt 
facilitates trapping at a multitude of weak, high-order mean motion resonances which, as 
Chiang et al. (2003) point out, is consistent with the detection of seven KBOs now known 
to librate at Neptune's 5:2 resonance. 

It is then natural to ask what mechanism might be responsible for stirring-up a broad 
swath of the Kuiper Belt, particularly since accretion models tell us that KBOs must have 
formed in a dynamically cold environment, i.e., the particles' initial e's and i's were < lO"^ 
(Stern 1996; Kenyon & Luu 1999). Note that this disturbance was probably not due to 
gravitational stirring by a number of long-gone protoplancts since, as Morbidelli et al. (2002) 
point out, a sizable fraction of such bodies would still persist in the Kuiper Belt and would 
likely have been discovered by now. 

Note that this stirring mechanism must also have had a large reach since it must afflict 
KBOs across the entire Main Belt, at least out to Neptune's 2:1 resonance. One mechanism 
that comes to mind is secular resonance sweeping, which is the only mechanism known to us 
that might stir eccentricities in the Belt up to e ~ 0.1 across its entire width (Nagasawa & Ida 
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2000). However this e-excitation is coherent in the sense that neighboring particles will have 
similar eccentricities. So it is unclear whether secular resonance sweeping of the Main Belt, 
which would then followed by sweeping mean-motion resonances due to Neptune's migration, 
will result in the range of eccentricities that is seen in Fig. 3. Secular resonance sweeping 
is a consequence of the dispersal of the solar nebula gas; the removal of that gas alters the 
giant planets' precession rates which in turn shifts the location of secular resonances (Ward 
1981). The magnitude of the disturbance caused by a sweeping secular resonance depends 
on the timescale over which the nebula is depleted; a longer depletion timescales Tdep results 
in larger eccentricity-pumping. The simulations of nebula dispersal by Nagasawa & Ida 
(2000) show that a disturbance of e ~ 0.1 across much of the Kuiper Belt requires a nebula 
depletion timescale of 10^ < Tdep ^ 10^ years (but see also Hahn & Ward 2002). 



13. Discussion 

One of the goals of this study is to determine how the adoption of a particular Kuiper 
Belt model might affect an assessment of the Belt's total population and mass. Note that 
some models of the Belt assume that the KBOs are distributed according to a primordial 
surface density distribution that might vary with distance at a{r) oc or so (e.g., Jewitt & 
Luu 1995; Trujillo et al. 2001), while the KBOs in another model are essentially equidistant 
(Bernstein et al. 2004). However Fig. 8 shows that a realistic Kuiper Belt has been eroded 
from the inside-out by the giant planets' gravitational perturbations, which suggests that 
the earlier models might not apply. However it turns out that an estimate of the total KBO 
population does not depend strongly upon a particular model's radial variation. As Section 
11 shows, all three models yield population estimates that are within a factor of ~ 2 of 
each other. This is because the observable KBOs really do inhabit a relatively narrow Belt 
centered on r ~ 45 AU having a radial half-width that is only Ar ~ 4 AU (see Fig. 8), 
so the the assumption of equidistant KBOs (e.g., Bernstein et al. 2004) appears to be good 
enough. 

It should also be noted that the magnitude interval over which a model Kuiper Belt 

can be compared to the observed Belt is given by the brightness of those KBOs having 
reliable orbits, and this sample is presently dominated by bodies having a relatively limited 
magnitude range of only 21 < m^j < 24. Further testing of this model, as well as the 
development of alternative models of the Belt, would be greatly facilitated if they could 
be compared to a larger sample of multi-opposition KBOs having reliable orbits and also 
exhibiting a broader range of apparent magnitudes and sizes. This larger KBO sample 
would be very useful in many ways. For example, it could be used to test the possibility 
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that the various Kuiper Belt subpopulations do exhibit variations in their size distributions 
(e.g., Section 6). This larger sample might also permit a better understanding of certain 
rare and unusual KBOs, such as those that inhabit the Extended Scattered Disk (Section 
10). A deeper understanding of the Kuiper Belt, and what the Belt tells us about the early 
evolution of the outer Solar System, would be facilitated by deeper KBO surveys over larger 
portions of the sky in a systematic way that leads to efficient KBO recoveries and reliable 
orbit determinations. 



13.1. summciry of findings 

• Accretion models have shown that Kuiper Belt Objects must have formed in a dynam- 
ically cold environment where the initial KBO seeds had nearly circular and coplanar 
orbits with eccentricities and inclinations < 0.001 (Stern 1996; Kenyon & Luu 1999). 
Simulations of Neptune's outwards migration into a dynamically cold Kuiper Belt, 
described in Section 3, show that the survivors in the Main Belt still maintain low 
eccentricities and inchnations. However this conflicts with the Main Belt's observed 
e's and sini's of ~ 0.1. This discrepancy suggests that some other process has also 
stirred-up the Kuiper Belt. This stirring event could have occurred prior to or after 
the onset of planet migration. 

• The existence of several KBOs hbrating at Neptune's 5:2 resonance suggest that this 
stirring event occurred prior to the onset of planet migration. Simulations by Chiang 
et al. (2003) have shown that if Neptune migrates into a stirred-up Kuiper Belt having 
eccentricities of e ~ 0.1, then trapping at Neptune's higher-order resonances, such 
as the 5:2, becomes more efficient. This result is confirmed by a higher-resolution 
study of this phenomena described in Section 4, which reveals that additional trapping 
also occurs at a number of exotic mean motion resonances like the 11:6, 13:7, 13:6, 
9:4, 7:3, 12:5, 8:3, 3:1, 7:2, and the 4:1; such resonances are not populated when 
Neptune migrates into a dynamically cold disk. Not surprisingly, Neptune's migration 
into a previously stirred-up Kuiper Belt also accounts for the eccentricities of e ~ 0.1 
observed in the Main Belt. 

• However the planet-migration scenario investigated here does not account for the ob- 
served KBOs having inchnations above i ~ 15° (Section 5), which is the main deficiency 
of this model. This is a serious discrepancy since half of all KBOs have inclinations 
i > 15° according to the debiased inchnation distribution reported by Brown (2001). 

• Neptune's migration into a stirred-up Kuiper Belt traps particles in eccentric orbits 
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at a number of resonances beyond a — 50 AU, the most prominent of these being the 
5:2 and the 3:1. Many of these distant particles that are trapped at semimajor axes 
a > 50 AU also have perihelia 30 < g < 40 AU, which is the domain conventionally 
known as the Scattered Disk. However Section 4 shows that only about 10% of the 
simulated particles that inhabit the so-called Scattered Disk or the Extended Scattered 
Disk (such as the gray zone in Fig. 4 where 50 < a < 80 AU and e > 0.25) are truly 
scattered particles. The vast majority of these particles never had a close encounter 
with Neptune; rather, they were placed in these wide, eccentric orbits by Neptune's 
sweeping mean-motion resonances. Note that the origin of these bodies as being due 
to resonant trapping is very distinct from the scattering scenario originally suggested 
by Duncan & Levison (1997). 

• Of the 10^ particles simulated here, only one managed to persist over the age of the 
Solar System in the Extended Scattered Disk, which is loosely defined as scattered 
orbits having periheha ? > 40 AU. This particle's orbit is qualitatively similar to 2000 
CRio5 which has a perihelion of g = 44 AU. However our simulations did not produce 
any extreme members of the Extended Scattered Disk that are similar to 2003 VB12 
(Sedna), which has a perihelion of g = 76 AU. 

• The output of the Nbody model is coupled to a Monte Carlo model that assigns radii R 
to the simulated particles according to a power-law type cumulative size distribution 
that varies as N{R) oc i?"'^. Magnitudes are computed for the simulation's particles, 
which then allows us to directly compare the simulated Belt to the observed Belt in a 
manner that accounts for telescopic selection effects. Section 6 compares the observed 
abundance of 2:1 objects to known Main Belt objects, and it is shown that the observed 
2:1 population is underdense by a factor of 20 relative to model predictions. Similarly, 
the observed 3:2 population is also depleted relative to model expectations. Another 
curious feature of the 3:2 is its lower than expected abundance (relative to the Main 
Belt KBOs) of fainter KBOs having magnitudes ttir > 23.5. Section 6 shows that this 
dearth of fainter KBOs at the 3:2 can be interpreted as a dearth of small bodies, which 
implies that the 3:2 population has a Q ~ 2.7 size distribution that is substantially 
shallower than the canonical Q — 4.4 power-law that holds for the larger members of 
the Main Belt. 

• The simulated Centaurs are quite sparse owing to their short dynamical lifetimes; only 
seven Centaurs were detected in the simulation during its final 2 Gyrs. Interestingly, 
all seven originated at/near Neptune's mean motion resonances in the Kuiper Belt 
(Section 8). The model puts an upper limit of Nc ^ 80 Centaurs having radii larger 
than R = 50 km, assuming they have a material density p = 1 gm/cm^ and an albedo 
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p — 0.04. It should be noted that all of the simulated Centaurs inhabit heliocentric 
distances of r > 22 AU, while the three Centaurs reported in Sheppard et al. (2000) 
were detected at r < 19 AU. If the simulated Centaurs are representative of reality, 
then this discrepancy in their heliocentric distances may indicate that Centaurs can 
break up and spawn additional Centaurs (e.g., Pittich & Rickman 1994) after evolving 
inwards from the Kuiper Belt. 

• This model also estimates that there are at most Nc ~ 1100 Trojans larger than 
i? = 50 km having a total mass of Mc ~ 5 x 10^^ orbiting at Neptune's triangular 
Lagrange points, assuming the usual p = 1 gm/cm^ and p = 0.04. 

• The absence of any distant KBOs having low eccentricities at a > 50 AU places tight 
upper limits on the abundance of any KBOs that might inhabit a hypothetical Outer 
Belt. Several upper limits arc inferred from this null result: (i) the primordial density 
of Outer Belt objects beyond 50 AU is smaller than the primordial Main Belt density 
by a factor / > 100, (m.) these distant KBOs are fainter than m* ~ 24.5 and thus 
have radii smaller than R ~ 80(p/0.04)~^/^ km, [in.) the cumulative size distribution 
of Outer Belt objects is steep, having a power-law index of Q > 6.0, or [iv.) some 
combination of the above. 

• The luminosity function of the Nbody/Monte Carlo model is fitted to the KBO's ob- 
served luminosity function, which then yields an estimate of the Belt's total population 
oi N 1.7 X 10^ KBOs larger than R ~ 50(p/0.04)~^/^ km having a total mass of 
Mtotai ~ 0.08(p/l gm/cm3)(p/0.04)--'^/2 ^^^g population and mass of the Belt's 
various subclasses (e.g.. Centaurs, Neptune Trojans, 3:2 and 2:1 populations, the Main 
Belt, and the Scattered Disk) are also assessed in Section 11 and hsted in Table 1. 
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Consider the orbital evolution of a particle trapped by a migrating planet's j + k : j 
mean-motion resonance. The planet's disturbing function contains the the resonant term 
Rres = Rres{c(, 6, (f)jk) whcrc a IS the planet/particle scmimajor axis ratio for a particle having 
an eccentricity e and a resonance angle having the form 



where A and Xp are the particle and planet's mean longitude while 00 is the particle's longitude 
of periapse (Murray & Dermott 1999). The relevant Lagrange planetary equations are 



A. Appendix A 
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where e = A — nt is the particle's mean longitude at epoch, n is the particle's mean motion, 
t is time, and coplanar orbits are assumed. According to the averaging principal, the non- 
resonant terms in the planet's disturbing function are of high frequency and average out 
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during the particle's libration period, and the planet's disturbing function is simply R — Rres- 
The derivatives of R are 

de d(j)jk duj d(j)jk 

With these, the planetary equations (A2) can be combined into a single differential equation, 

' VI -e2 - — . (A4) 



dt 2e \j + k J a dt 



It will be convenient to replace e with the variable E = \/\ — e?'\ since ede = —EdE, the 
above equation can then be recast as 

da 2dE , 

where 7 = j/{j + k). This differential equation is now easily integrated and yields Ina = 
—2 ln(7 — E') + In 5 where the integration constant B can also be expressed as 

B^a( - (A6) 

V J + kJ 

{c.f., Hahn & Ward 1995). The earliest derivation of this integral of the motion known to 
us is given in Brouwer (1963). 

Brouwer originally used this integral to consider the motion of a particle being perturbed 
by a planet in a static orbit. However Yu & Tremaine (1999) recognized that this integral 
is preserved even when the planet is migrating, which means that Eqn. (A6) can be used to 
predict the particle's eccentricity as its orbit is expanded by the planet's resonance: 

If, prior to capture, the particle is in a circular orbit with an initial semimajor axis ao, then 
B = ao[k/{j + k)]^ and 

e\a)^l-(^-±^^J , (A8) 

and the particle's eccentricity grows as the particle's orbit expands. An approximate form 
of this expression is also derived in Malhotra (1993b). Thus if Neptune's orbital expansion 
is indeed responsible for the KBOs seen at Neptune's 3:2 resonance at a = 39.4 AU having 
eccentricities as high as e ~ 0.32 (see Fig. 2), Eqn. (A8) indicates they must first have been 
orbiting at ao = 28.0 AU at the time of capture, and that Neptune was initially orbiting at 
ajv = OoT^''^ — 21.4 AU and hence migrated a distance Ajv = 8.7 AU. 



-34- 



B. Appendix B 

The following will show that the ratio of the observed abundance of any two dynamical 
classes of KBOs is approximately equal to the ratio of the absolute abundances of the much 
larger unseen populations. Begin by letting N°^'^{m) be the number of observed KBOs that 
inhabit some dynamical class x that are brighter than apparent magnitude m, where x might 
represent, say, the 3:2 population. Also assume that all KBOs have the same heliocentric 
distance r, which is a common assumption (e.g., Irwin et al. 1995; Bernstein et al. 2004) 
that simplifies this analysis considerably, and is actually not a bad assumption since Fig. 8 
shows that most of the observed multi-opposition KBOs considered here do indeed inhabit a 
rather narrow belt, with more that half of these KBOs found at heliocentric distances within 
3 AU of r = 42 AU. We will also assume that the Kuiper Belt has azimuthal symmetry, 
which is justified further in Appendix C. We shall also assume that all the KBO astronomers 
are observing this Belt largely along the ecliptic, which also is not a bad assumption, since 
72% of the multi-opposition KBOs studied here have latitudes of < 3°, which is the 
typical inclination of the Belt's low-i component (Brown 2001). This indicates that the 
astronomers who discovered most of the KBOs in this sample (represented by the red dots 
in Fig. 3) were observing largely along the Belt's midplane. So, although these assumptions 
are not rigorously correct in detail, they are good enough to allow us to assess qualitatively 
the relative abundances of the various KBO populations. 

Let S-r(m) represent the sky-plane number density of class-x KBOs {i.e., the cumulative 
luminosity function of class-x KBOs), while S^(m) = dLx/dm be the differential luminosity 
function for class x. Then number of observed KBOs brighter than magnitude m is (Bernstein 
et al. 2004) 

iVf ^(m) = j E'^{m)n{m)ri^{m)dm, (Bl) 

where Q{m) is the total sohd angle that has been surveyed by all KBO astronomers to a 
limiting magnitude m, and r]x{m) is the efficiency at which a KBO of class x and magnitude 
m and is detected. Thus the number of observed class x KBOs of having magnitudes in the 
interval m ± Am/2 is 

dN?^\m) = —^Am = J:'(m)n(m)r]Jm)Am. (B2) 
dm 

If we were to compare class x KBOs to, say, class y KBOs, then the ratio of their observed 
abundances is 

- d^) " eJR' 

where it is assumed that the KBO detection efficiency is insensitive to dynamical class, i.e., 
rix{m) ~ Vyi"^)- Note that this ratio depends only on the KBOs' differential luminosity 
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functions S'(m), and that it is insensitive to their discovery details such as the total solid 
angle Q(m) that the KBO astronomers have surveyed to depth m. 

The cumulative luminosity function for class x can be written as E^(m) = Nx{m)/fltotai 
where ^totai is the Kuiper Belt's total solid angle, and Nx{m) is the cumulative magnitude 
distribution for class x, ie, the total number of KBOs in class x that are brighter than 
magnitude m. Thus 

E' r 1 = = ^—^ (ba) 

"^""^ ^totai \dmj ntotai \dRj \dmj ^ ' 



since a body's magnitude m is a function of its radius R via Eqn. (6), for which 

(B5) 



dR Rfnin \ Rfni 

is the differential size distribution for class x which has a total of N^.totai bodies having radii 
in the interval Rmin < R < Rmax (see Eqn. 4). And since dR/dm = —R/5 (see Eqn. 6), the 
differential luminosity function can be written 

SL(m(i?)) = f-^V^ (B6) 



Thus if any two KBO classes x and y have the same size distribution Q, then the ratio of 
their observed abundances is simply 

Y] [m) N.J. total 

2:'y{m) Ny^tOtal 

In other words, the ratio of the observed abundance of any two classes of KBOs is approx- 
imately equal to the ratio of their intrinsic abundances, provided all bodies have the same 
size distribution Q. But if populations x and y have different size distributions, then rx:y{m) 
will vary with m (see Section 6). 



C. Appendix C 

The Monte Carlo model of Section 6 replicates each Nbody survivor seen in Fig. 3 10^ 
times by randomizing the particles' mean anomalies, which results in a model Kuiper Belt 
having azimuthal symmetry. Thus the model implicitly assumes that the Kuiper Belt has 
this same symmetry. However one might question this assumption, since it is well known that 
a planet's resonant perturbations can rearrange the longitudes of a small body population. 
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For instance, particles at Neptune's 3:2 tend to approach perihelion at longitudes that are 
~ 90° away from the planet (Malhotra 1995). This is illustrated in the upper left corner of 
Fig. 16, which shows the ecliptic coordinates of all simulated particles that are in or very near 
Neptune's 3:2 resonance. This figure shows that Neptune tends to arrange the 3:2 bodies 
preferentially away from the Sun-Neptune line, causing regions that lead/trail Neptune by 
±90° to be more densely populated, termed 'sweet spots' in Chiang & Jordan (2002). Figure 
16 shows that resonant shepherding also occurs at Neptune's 2:1, again similar to that seen 
in Chiang & Jordan (2002). Red dots indicate the positions of the multi-opposition KBOs 
considered here, which shows that KBOs in or near the 3:2 resonance tend to be discovered 
at longitudes that are roughly ±90° away from Neptune, as expected. 

It is thus possible that Neptune's rearrangement of the Belt might skew our estimate 
of the resonant populations. If, for example, KBO astronomers were preferentially sur- 
veying the Belt along longitudes that are ±90° away from Neptune, they would detect a 
rather high sky-plane number density of KBOs in the 3:2. If we then assumed that this 
3:2 number density were uniform about the entire ecliptic, we would overestimate the total 
3:2 population. Similarly, if astronomers systematically observed the Belt towards Neptune 
or 180° away, we would underestimate the 3:2's average sky-plane density and undercount 
its total population. But if astronomers surveyed all longitudes with equal frequency, then 
these competing effects — due to Neptune pushing KBOs away from certain longitudes and 
towards other longitudes — should wash out, resulting in an estimate of the 3:2 population 
that is approximately reliable. 

The upper right part of Fig. 16 shows the ecliptic coordinates of the observed multi- 
opposition Main Belt KBOs. Since the Main Belt population is likely azimuthally symmetric, 
these red dots should be a good indicator of where astronomers are looking for KBOs. It 
is quite clear from this figure that these astronomers' lines-of-sight are not distributed 
uniformly about the ecliptic. For instance, these KBO hunters tend to avoid the galactic 
plane, which passes through the ecliptic along the dashed line. The MB figure also shows 
that astronomers have avoided a narrow portion of the 3:2's 'sweet spot' that leads Neptune 
by about 90°. But it is also quite clear that these same astronomers are not systematically 
staring at the 3:2's sweet spots, nor are they systematically avoiding them. Consequently, our 
simple visual inspection of this KBO sample indicates that our estimate of the 3:2 population 
will not be significantly baised towards under or overcounting the 3:2 population, and that 
the 3:2 abundance reported in Section 11 is indeed representative. A similar conclusion is 
also drawn for the much sparser 2:1 population. 

Of course, this treatment succeeds only if the Monte Carlo method of Section 6 does 
not alter the particles' radial distribution as they have their mean anomalies randomized by 
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the replication process. However it is straighforward to show that this is indeed the case. 
For instance, randomizing the mean anomahes of the 3:2 population seen Fig. 16 does not 
change their radial distribution in any significant way. 
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Table 1. Abundance and mass of KBOs having radii 50 < i? < 1000 km assuming density 

p—1 gm/cm^ and albedo p — 0.04 



dynamical class 


site (AU) 


population 


mass (M( 


b) 


Centaurs 


a < ajsi 


< 80 


< 4 X 10- 


-5 


Trojans 


a = ojv 


< 1.1 X 10^ 


< 5 X 10 


-4 


3:2 


a = 39.5 


2.7 X 10^ 


3.1 X 10- 


3 


Main Belt 


40.1 <a< 47.2 


1.3 X 10^ 


0.059 




2:1 


a = 47.8 


5.3 X 10^ 


2.4 X 10" 


-3 


Scattered Disk 


50 < a < 150 & 28 < g < 40 


2.5 X 10^ 


0.011 
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Fig. 1. — The planets' semimajor axes versus time t, as well as a few of Neptune's mean 
motion resonances that are effective at trapping particles. 
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Fig. 2. — A simulation of Neptune's migration into a dynamically cold Kuiper Belt. The 
model Kuiper Belt is initially composed of 10^ massless particles having semimajor axes 
randomly distributed over 20 < a < 80 AU with a surface number density that varies as 
a~^. The particles initial eccentricities are Rayleigh distributed about a mean (e) = 0.001, 
and their initial inclinations (which are measured with respect to the system's invariable 
plane) are similarly distributed over (sini) = (e)/2, while their other angular orbit elements 
are uniformly distributed over 27r. The four giant planets' orbits migrate according to the 
prescription described in Section 2, and the black dots shows the simulated Kuiper Belt 
endstate — the particles' e's and i's versus a at time t = 5 x 10* years. The red dots are 
the ecliptic orbit elements of KBOs reported by the Minor Planet Center as having been 
observed for at least two oppositions. The yellow dots indicate Neptune's orbit, and the 
vertical dashes show the locations of Neptune's various j + k : j mean-motion resonances. 
Orbits having periheha q — 30 and 40 AU are also indicate by the curves. 




<-*lk kk kki-A **-k-k kk ■■■■ 

iTJ^ ro inr^ cm tnr^ in d) in 




30 40 50 50 70 80 

seiTfinnajor axis a (AU) 

Fig. 3. — Neptune migrates into a dynamically hot disk of = 10^ particles having initial 
e's Rayleigh distributed about (e) = 0.1 and inclinations similarly distributed about (sini) = 
(e)/2, while all other initial orbits are distributed as in Fig. 2. This system is evolved for 
t = 4.5 X 10^ years. 
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Fig. 4. — Eccentricity e versus semimajor axis a at time t = 4.5 x 10^ years for the model 
population having initial e ~ 0.1 shown in Fig. 3. The vertical dashes indicate the mean 
motion resonances that are occupied by trapped particles with perihelia g < 40 AU or 
particles inhabiting resonances in the gray region with libration amplitudes \4>jk\ < 90°. 
Crosses indicate scattered particles. 
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Fig. 5. — The ecliptic inclination distributions for the simulated and observed bodies in Fig. 
3 that have perihelia g < 42 AU and latitudes < 1°. Error bars give the Poisson counting 
uncertainty, and the vertical scale is arbitrary. 
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Fig. 6. — The abundances of the Belt's various dynamical populations are plotted over time t, 
with all populations normahzed to the Main Belt (MB) population at time i = 4.5 x 10^ years, 
and counting only those particles having perihelia g < 45 AU. The resonant populations are 

counts of all particles having semimajor axes ±0.6 AU within exact resonance, while the MB 
curve gives the number of particles having semimajor axes 40.1 < a < 47.2 AU, excepting 
bodies in/near the 5:3 and 7:4 resonances. However the Scattered Disk (SD) curve is the 
number of particles having a > 49 AU and perihelia q < AO AU that are not members of any 
of the indicated resonances (e.g., the 5:2 and 3:1), while Centaurs (C) are particles having 
a < Neptune's semimajor axis. 
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Fig. 7. — The dark-gray curve in the upper plot is the number of KBOs that have been 
observed (for two or more oppositions) orbiting in/near the 2:1 resonance having apparent 
magnitudes m/?± Am relative to the observed Main Belt (MB) population at 40.1 < a < 47.2 
AU, while the lower dark-gray curve is the observed 3:2/MB ratio. Likewise, the light-gray 
curves are the 2:1/MB and 3:2/MB ratios for the simulation of Fig. 3, where these bodies' 
sizes and magnitudes are assigned using the Monte Carlo method of Section 6 assuming 
Q — 4.4 and Rmin — 20 km. All of these curves are smoothed over a magnitude-window 
having a half-width Am = 0.5, and the vertical half-width of the gray zones are one standard 
deviation assuming Poisson counting uncertainties. The observed curves end at rriR ~ 24.5, 
which is the magnitude of the faintest multi-opposition KBO. The dashed curve is the 
simulated 3:2/MB ratio when the 3:2 population has a shallow size distribution with Q = 2.7 
and Rmin = 4.3 km, while the MB bodies have Q = 4.4 and Rmin = 20 km. 
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Fig. 8. — Dots show the KBO surface number density o"(r), normahzed to peak at unity, that 
is inferred from the KBOs radial distribution function /(r) oc ra{r) reported by Trujillo & 
Brown (2001). The t — and 4.5 Gyr curves are the simulated Belt's initial and final surface 
density. The 'truncated Belt' curve is the simulation's final surface density assuming that 
the Belt is truncated at a = 45 AU, and that the 3:2 population also does not contribute 
to a. Note that the surface density of the inner half of this Belt increases as a very steep 
function of distance as a{r) oc r^^^ for r < 45 AU. 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




Fig. 9. — The gray curve is tiie expected ratio of bodies in tiie Outer Belt (OB, beyond 
a > 50 AU) relative to the Main Belt (MB, 40.1 < a < 47.2 AU) as a function of magnitude 
niR assuming the size distribution has Q — 4.4 and Rmin — 20 km. The solid curve is the 
OB/MB ratio assuming the OB instead has a steeper Q — 6.0 size distribution while the 
MB has Q — 4.4, with Rmin — 20 km for both populations. 




Fig. 10. — Open circles are the orbits of the simulation's seven Centaurs at time t = 10^ 
years, while black dots indicate their subsequent motions. Gray dots are the simulation's 
endstate at time t = 4.5 x 10^ years, red dots indicate multi-opposition KBOs and Centaurs, 
and yellow dots are the giant planets' final orbits. Neptune's mean-motion resonances are 
indicated, and the curve is the threshold for Neptune-crossing orbits. 
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Fig. 11. — The black dots are the longitudes of five simulated trojans versus their semimajor 
axes, all relative to Neptune's, and sampled every 10^ years during the entire simulation. 
Neptune lies at the +. Small squares show the final positions of these trojans, indicating 
that two trojans he at the leading Lagrange point, with three trailing. Red dots indicate 
three nearby 'field' KBOs as well as Neptune's two known Trojans, 2001 QR322 and 2004 
UPio. 
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Fig. 12. — Eccentricities e are plotted versus a on a logarithmic axis for the particles of Fig. 
3. Red dots are the KBOs observed over multiple oppositions, and the curves indicate the 
g = 40 AU threshold as well as Neptune-crossing orbits. Two members of the so-called 
extended Scattered Disk, 2000 CR105 and 2003 VB12, are also indicated. Crosses indicate 
particles resonantly trapped at the 3:1, 7:2, and the 4:1. The large black dot indicates the 
only simulated particle that was scattered into a high-perihelia orbit {q = 43.5 AU) that 
resembles 2000 CR105. 
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Fig. 13. — The orbital history of the only simulated scattered object that inhabits the 
extended Scattered Disk (e.g., the large black dot in Fig. 12). The particle's a, e, g , i, and 
Tisserand parameter T — aN / a-\-'2^ {a/aN){}- — e^) cosi are first plotted versus logarithmic 
time (left half of the Figure) and then linearly (right half). The dotted curves are Neptune's 
semimajor axis ajv and a few of its mean-motion resonances. 




Fig. 14. — The dark gray curve is the observed ratio of Scattered Disk (SD) objects in orbits 
having 50 < a < 150 AU and perihelia 28 < 5 < 40 AU, relative to Main Belt (MB) objects, 
versus apparent magnitude tur. The light gray curve is the simulated ratio. 
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Fig. 15. — The smooth black curve is the Kuiper Belt's observed ecliptic luminosity function 
E(mij), obtained by integrating the differential luminosity function dTi{mfj) / dniR reported 
in Bernstein et al. (2004); note that it breaks to a shallower slope at magnitudes rribreak — 24. 
Dots give the simulation's ecliptic luminosity function for particles in a truncated Kuiper 
Belt (see Section 11) having latitudes within 0.5° of the ecliptic. Error bars are for Poisson 
counting statistics. 
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Fig. 16. — Black dots indicate the positions of all simulated particles relative to Neptune 
recorded at all times later than t = 3 Gyrs. Red dots indicate the ecliptic coordinates of 
the observed multi-opposition KBOs, shown for July 1, 2000, which is the date by which 
half of the sample considered here had been discovered. The blue dot is Neptune's position 
on this date, and the dashed line shows where the galactic plane penetrates the ecliptic x-y 
plane. The particles and KBOs are sorted by their dynamical membership: bodies in or very 
near the 3:2 resonance (upper left). Main Belt (MB, upper right), 2:1 (lower left), and the 
so-called Scattered Disk Objects (SD, lower right) having perihelia 30 < g < 40 AU and 
a > 48.4 AU. 



